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ABSTRACT 


We  continue  exploring  methodologies  to  improve  earthquake-explosion  discrimination  using  regional 
amplitude  ratios  such  as  P/S  in  a  variety  of  frequency  bands.  Empirically,  we  demonstrate  that  such  ratios 
separate  explosions  from  earthquakes,  using  closely  located  pairs  of  earthquakes  and  explosions  recorded  on 
common,  publicly  available  stations  at  test  sites  around  the  world  (e.g.,  Nevada,  Novaya  Zemlya, 
Semipalatinsk,  Lop  Nor,  India,  Pakistan,  and  North  Korea).  We  are  also  examining  if  there  is  any  relationship 
between  the  observed  P/S  and  the  point  source  variability  revealed  by  longer  period  full  waveform  modeling 
(e.g..  Ford  et  ah,  2008).  For  example,  regional  waveform  modeling  shows  strong  tectonic  release  from  the  May 
1998  India  test,  in  contrast  with  very  little  tectonic  release  in  the  October  2006  North  Korea  test,  but  the  P/S 
discrimination  behavior  appears  similar  in  both  events  using  the  limited  regional  data  available. 

While  regional  amplitude  ratios  such  as  P/S  can  separate  events  in  close  proximity,  it  is  also  empirically  well 
known  that  path  effects  can  greatly  distort  observed  amplitudes  and  make  earthquakes  appear  very  explosion¬ 
like.  Previously  we  have  shown  that  the  Magnitude  Distance  Amplitude  Correction  (MDAC)  technique 
(Walter  and  Taylor,  2001)  can  account  for  simple  1-D  attenuation  and  geometrical  spreading  corrections,  as 
well  as  magnitude  and  site  effects.  However,  in  some  regions,  1-D  path  corrections  are  a  poor  approximation, 
and  we  need  to  develop  2-D  path  corrections.  Here  we  demonstrate  a  new  2-D  attenuation  tomography 
technique  using  the  MDAC  earthquake  source  model  applied  to  a  set  of  events  and  stations  in  both  the  Middle 
East  and  the  Yellow  Sea  Korean  Peninsula  regions.  We  believe  this  new  2-D  MDAC  tomography  has  the 
potential  to  greatly  improve  earthquake-explosion  discrimination,  particularly  in  tectonically  complex  regions 
such  as  the  Middle  East. 
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OBJECTIVES 


Monitoring  the  world  for  potential  nuclear  explosions  requires  characterizing  seismic  events  and  discriminating 
between  natural  and  man-made  seismic  events,  such  as  earthquakes  and  mining  activities,  and  nuclear  weapons 
testing.  We  continue  developing,  testing,  and  refining  size-,  distance-,  and  location-based  regional  seismic 
amplitude  corrections  to  facilitate  the  comparison  of  all  events  that  are  recorded  at  a  particular  seismic  station. 
These  corrections,  calibrated  for  each  station,  reduce  amplitude  measurement  scatter  and  improve  discrimination 
performance.  We  test  the  methods  on  well-known  (ground  truth)  datasets  in  the  U.S.  and  then  apply  them  to  the 
uncalibrated  stations  in  Eurasia,  Africa,  and  other  regions  of  interest  to  improve  underground  nuclear  test 
monitoring  capability. 


RESEARCH  ACCOMPLISHED 


As  part  of  the  overall  National  Nuclear  Security  Administration  Ground-Based  Nuclear  Explosion  Monitoring 
Research  and  Development  (GNEMRD)  program,  we  continue  to  pursue  a  comprehensive  research  effort  to 
improve  our  capabilities  to  seismically  characterize  and  discriminate  underground  nuclear  tests  from  other  natural 
and  man-made  sources  of  seismicity.  To  reduce  the  monitoring  magnitude  threshold,  we  make  use  of  local  and 
regional  body  and  surface  wave  data  to  calibrate  each  seismic  station.  Our  goals  are  to  reduce  the  variance  and 
improve  the  separation  between  explosions  and  other  event  populations  by  accounting  for  the  effects  of  propagation 
and  differential  source  size  and  by  optimizing  the  types  and  combinations  of  amplitude  measurements  used. 


Empirical  P/S  Observations 

Empirical  observations  have  shown  clear  differences  between  earthquake  and  explosion  P/S  amplitude  ratios  that 
can  be  exploited  for  nuclear  explosion  monitoring.  These  differences  became  apparent  with  the  advent  of 
widespread  digital  seismometry  in  the  late  1980s  and  early  1990s,  when  a  large  number  of  different  researchers 
found  that  at  high  frequencies,  greater  than  around  2-4  Hz,  the  regional  P/S  amplitude  ratios  consistently 
discriminated  explosions  from  earthquakes  (e.g.,  Dysart  and  Pulli,  1987;  Baumgardt  and  Young,  1990;  Kim  et  ah, 
1993;  Walter  et  ah,  1995;  Taylor,  1996;  Hartse  et  al.,  1997;  Battone  et  ah,  2002).  For  example,  in  Figure  1  we 
compare  an  earthquake  and  an  explosion  at  each  of  four  major  test  sites  (rows),  bandpass  filtered  in  four  different 
frequency  bands  (columns).  Here  we  use  closely  spaced  pairs  of  earthquakes  and  explosions  recorded  at  the  same 
station  to  minimize  any  path  and  site  effects,  and  so  the  observed  differences  come  from  effects  in  the  source 
region.  There  is  a  relative  lack  of  S-waves  in  the  explosions  as  the  frequency  increases.  Note  also  there  are 
interesting  differences  between  the  test  sites,  indicating  that  emplacement  conditions  (depth,  media  properties, 
geologic  and  tectonic  conditions,  topography,  etc.)  play  a  role  in  relative  P-  and  S-wave  generation. 

Examination  of  P/S  ratios  at  high  frequencies  has  now  been  done  all  over  the  world.  As  we  show  in  Figure  2,  the 
characteristic  pattern — that  explosions  have  larger  relative  P-wave  amplitudes  than  S-wave  amplitudes — appears  to 
hold  everywhere,  provided  that  we  filter  the  data  at  a  high  enough  frequency  and  that  we  compare  events  with  very 
similar  paths.  One  example  in  Figure  2  is  the  announced  nuclear  test  by  North  Korea  on  October  9,  2006.  In  this 
case,  the  nearest  earthquakes  for  comparison  are  around  100  km  away,  but  we  can  bracket  the  test  with  earthquakes 
and  clearly  see  that  it  is  deficient  in  high  frequency  S-wave  energy.  This  example  shows  the  importance  of 
accounting  for  path,  as  the  S-wave  energy  in  this  case  is  partitioned  into  both  Sn  and  Eg  phases.  The  path  of  the 
recorded  seismic  energy  crosses  the  Sea  of  Japan,  a  region  of  oceanic  cmst  where  Eg  does  not  efficiently  propagate. 
On  can  see  that  this  earthquake,  with  a  greater  proportion  of  its  path  in  the  Sea  of  Japan,  has  stronger  Sn  than  the 
other.  It  is  important  to  correct  for  path  effects  and  to  consider  P/S  ratios  for  all  the  observable  regional  phases  and 
avoid  cases  where  one  of  the  earthquake  S-wave  phases  is  attenuated  to  below  the  level  of  the  noise  and  it  appears 
explosion-like  simply  due  to  the  path  effects. 

When  similar-sized  earthquakes  and  explosions  are  nearly  co-located,  we  can  understand  the  observed  seismic 
contrasts,  such  as  the  relative  P-to-S  wave  excitation,  in  terms  of  depth,  material  property,  focal  mechanism,  and 
source  time  function  differences.  However,  it  is  well  known  that  path  propagation  effects  (e.g.,  attenuation, 
blockage)  and  source  scaling  effects  (e.g.,  comer  frequency  scaling  with  magnitude)  can  make  earthquakes  look  like 
explosions  and  vice  versa.  We  developed  the  MDAC  technique  (Walter  and  Taylor,  2001)  to  account  for  these 
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Fignre  1.  An  earthqnake  and  explosion  pair  at  fonr  different  test  sites:  Nevada,  U.S.A.;  Novaya  Zemiya, 
Rnssia;  Semipaletinsk  former  U.S.S.R.  now  Kazakhstan;  and  Lop  Nor,  China.  The  seismograms 
are  bandpass  filtered  at  fonr  freqnency  hands  (colnmns)  illnstrating  significant  differences  in 
relative  P  and  S  wave  amplitndes,  with  explosions  have  relatively  less  S-waves  at  higher  freqnencies. 

effects  with  proper  calibration.  We  use  the  earthquakes  alone  to  determine  the  MDAC  parameters,  such  as 
geometrical  spreading,  frequency  dependent  Q,  and  the  average  apparent  stress.  After  calibration,  the  MDAC 
formulation  provides  expected  spectral  amplitudes  as  a  function  of  phase,  frequency,  magnitude  and  distance.  These 
can  then  be  subtracted  from  the  actual  observations.  For  earthquakes,  the  corrected  data  should  exhibit  a  close-to- 
zero  mean  and  no  significant  trends  with  magnitude  and  distance.  Explosions  should  have  significant  non-zero 
mean  residuals,  leading  to  improved  discrimination. 

As  shown  in  Figure  3,  we  applied  MDAC  to  correct  for  path  effects  for  earthquakes  all  around  the  Korean 
Peninsula,  Yellow  Sea,  and  northeastern  China  region  recorded  at  stations  MDJ  and  TJN,  which  are  publicly 
available  through  data  centers  such  as  the  Incorporated  Research  Institutions  for  Seismology  (IRIS)  and  the  Ocean 
Hemisphere  Project  Data  Management  Center  (OHP-DMC).  When  ratios  were  available  from  both  stations,  they 
were  averaged.  We  look  at  three  different  P/S  ratios  in  four  frequency  bands.  As  Figure  3  shows,  the  North  Korean 
test  appears  to  discriminate  progressively  better  at  higher  frequencies  for  all  three  types  of  P/S  ratios.  Due  to  event 
size,  attenuation,  background  noise,  and  the  need  to  be  past  the  crossover  distance  to  see  Pn  and  Sn,  not  all 
measurements  can  be  made  for  all  phases  and  frequency  bands.  There  are  167  earthquakes  for  1-2  Hz  Pg/Lg  and 
fewer  for  the  other  ratios.  In  general,  the  number  of  measurements  passing  the  signal-to-noise -ratio  (SNR) 
thresholds  (here  2  for  P-waves  and  1.3  for  S-waves)  diminishes  as  the  frequency  increases. 
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High-frequency  P/S  appears  to  work  everywhere 
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Fignre  2.  Fonr  examples  of  6-8  Hz  filtered  seismograms  of  explosions  compared  to  nearhy  earthqnakes 
ontside  the  major  testing  areas  show  the  characteristic  difference:  explosions  have  larger  P/S 
amplitnde  valnes. 

Practical  considerations,  such  as  whether  we  expect  a  new  event  in  a  particular  area  with  a  given  magnitude  to  be 
detectable  for  a  given  phase  and  frequency,  are  critical  for  monitoring.  Therefore  it  is  good  to  see  that  there  remains 
considerable  discrimination  capability  in  all  three  types  of  P/S  ratios  in  the  bands  below  6-8  Hz,  particularly  the  2-4 
and  4-6  Hz  bands.  As  we  have  noted  in  prior  MRR  papers,  it  is  possible  to  form  multivariate  combinations  of  ratios 
to  enhance  discrimination,  and  these  can  be  tuned  to  optimize  discrimination  based  on  the  bands  you  are  likely  to 
have  to  work  with. 

One  type  of  ratio  that  has  been  considered  from  a  practical  viewpoint  is  the  cross-spectral  ratio  (e.g.,  Taylor  et  al., 
2002):  comparing  a  higher  frequency  P-wave  amplitude  to  a  lower  frequency  S-wave  amplitude.  As  shown  in 
Figure  1  for  Lop  Nor  or  Figure  2  for  India,  some  explosions  do  not  have  any  visible  S-waves  at  high  frequency  and 
the  S-measurement  is  made  on  coda.  Such  explosion  coda  measures  can  work  so  long  as  one  is  sure  a  comparable- 
size  earthquake  would  have  produced  measurable  S-waves  allowing  discrimination.  In  a  cross-spectral  ratio,  the  S- 
wave  band  is  chosen  so  that  even  explosion  S-waves  are  detectable,  and  these  have  been  shown  to  work  well  in 
some  cases  such  as  Lope  Nor  (e.g.,  Taylor  et  al.,  2002).  A  complication  with  cross-spectral  ratios  occurs  if  the 
spectral  decay  at  high  frequencies  is  substantially  different  for  explosions  than  for  earthquakes:  it  can  start  to 
counteract  the  P/S  difference  and  lead  to  poor  discrimination.  In  Figure  4,  we  show  an  example  of  a  6-8  Hz  Pn  to 
2-4  Hz  Lg  ratio  for  Nevada  Test  Site  (NTS)  data  averaged  over  two  stations,  ELK  and  KNB.  One  can  see  that  a 
number  of  the  explosions  do  not  discriminate  at  all.  NTS  explosions  show  a  clear  material  property  dependence 
(e.g.,  Walter  et  al.,  1995)  with  explosions  in  weak,  high-gas-porosity  media  having  much  stronger  spectral  decay 
than  explosions  in  stronger,  low-gas-porosity  media  when  compared  at  the  same  station  using  the  same  path  and  site 
corrections,  as  shown  in  Figure  4. 
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Figure  3.  The  map  shows  earthquakes  (blue  circles)  aud  the  October  9,  2006,  North  Koreau  nuclear  test  (red 
star)  observed  at  seismic  stations  MDJ  and  TJN.  The  scatter  plots  show  the  MDAC  path  corrected 
P/S  ratios  at  each  station  (average  when  both  available)  for  three  different  P/S  ratios  in  four 
different  frequency  bands. 
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Figure  4.  Left:  example  cross-spectral  ratio  for  NTS  area  data  averaged  for  stations  ELK  and  KNB.  Right: 

NTS  explosions  show  strong  material  property  dependence  of  the  P  and  S  spectra.  We  separate  the 
nuclear  explosions  into  high-strength/low-gas-porosity  (black)  and  low-strength/high-gas-porosity 
(red)  populations,  following  the  criteria  in  Walter  et  al.  (1995)  of  log(pa^)  =  0.093GP  +  9.4. 
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In  Figure  5,  we  show  the  same-frequeney  P/S  ratios  in  different  frequeney  bands  for  the  ELK-KNB  station  averaged 
data  eolor  eoded  as  in  Figure  4  by  media  properties  from  Springer  et  al.  (2002).  At  lower  frequeneies  (0.5-4  Hz) 
there  is  some  evidenee  that  the  material  properties  affeet  diserimination  with  the  weak,  high  gas  porosity  events  do 
not  diseriminate  as  well.  At  higher  frequeneies  (6-16  Hz)  there  does  not  appear  to  be  a  signifieant  differenee  though 
the  dataset  starts  to  get  fairly  small.  Higher  frequeney  P/S  ratios  (6  Hz  and  above)  have  beeome  an  important 
diseriminant  when  they  are  available  beeause  of  this  apparent  laek  of  material  property  dependenee  and  good 
diserimination  performanee  in  all  regions  where  good  empirieal  data  is  available. 


Pg/Lg  at  ELK-KNB  path  corrected  with  MDAC 


Myy  M\fj  Myy 


Figure  5.  Pg/Lg  ratios  for  two-station  (ELK-KNB),  averaged,  MDAC-corrected,  NTS-region  data.  We 

separate  the  nnclear  explosions  into  high-strength/low-gas-porosity  (hlack)  and  low-strength/high- 
gas-porosity  (red)  popnlations,  as  in  Figure  4.  There  is  some  evidence  of  a  material  property  effect 
at  low  frequencies,  hut  not  at  frequencies  above  6  Hz.  Note  that  discrimination  performance  in 
general  continues  to  improve  as  frequency  increases,  at  least  to  12  Hz,  hut  the  number  of  events 
passing  the  SNR  decreases  significantly. 


Are  Long-Period  Complexities  Reflected  in  Short-Period  Discriminants? 

A  remarkable  feature  of  the  high  frequency  P/S  ratios  for  discriminating  underground  nuclear  tests  is  that  they  seem 
to  work  about  as  well  in  cases  where  the  long-period  signals  indicate  that  the  explosions  are  complex  as  in  cases 
where  the  long-period  signals  are  very  simple.  For  example,  if  we  contrast  the  May  1 1,  1998,  Indian  test  with  the 
October  9,  2006,  North  Korean  test,  both  show  good  high-frequency  P/S  discrimination  in  Figure  2.  However,  the 
Indian  test  shows  significant  Love  waves  and  reversed  Rayleigh  waves  (Walter  et  al,  unpublished  data;  Selby  et  al., 
unpublished  data),  indicating  a  large  amount  of  tectonic  release  and  significant  non-isotropic  components  of  the 
moment  tensor  source.  On  the  other  hand,  waveform  modeling  of  the  North  Korean  test  shows  almost  no  Love 
waves  and  a  nearly  pure  isotropic  moment  tensor  (e.g.,  Walter  et  al,  2007).  Despite  the  complexities  in  the  Indian 
test  source,  the  S-waves  disappear  at  frequencies  of  1  Hz  and  higher,  and  the  test  easily  discriminates  from  nearby 
earthquakes,  as  shown  in  Figure  2  and  by  Rogers  and  Walter  (2002). 

To  better  understand  any  connections  between  complexities  in  the  long  period  signal  and  the  high  frequency 
discriminants,  we  are  making  use  of  the  waveform  modeling  results  of  Ford  et  al.,  2008  (see  also  Ford  et  al.,  these 


707 


2008  Monitoring  Research  Review:  Gronnd-Based  Nnclear  Explosion  Monitoring  Technologies 


Proceedings).  For  example,  it  has  been  proposed  that  the  explosions  can  have  a  significant  compensated  linear 
vector  dipole  (CLVD)  component  as  part  of  their  interaction  with  the  free  surface  and  that  this  can  impact  surface 
wave  magnitude  and  perhaps  S-wave  amplitudes  as  well  (Patton  and  Taylor,  2008).  In  Figure  6,  we  compare  Pg/Lg 
ratios  in  two  frequency  bands  to  the  relative  CLVD  strength  from  waveform  modeling.  At  the  1-  Hz  band  there  may 
be  a  correlation  that  is  not  apparent  at  the  high  frequencies.  However,  we  are  still  assessing  just  how  well  resolved 
the  relative  CLVD  components  from  the  waveform  model  are  before  reaching  any  conclusions.  We  hope  to  have 
updated  results  to  present  at  the  meeting.  Finally,  we  note  that  there  is  clear  evidence  from  the  coda  that  the  absolute 
depth  of  burial  has  an  affect  on  Rg  excitation  (e.g.,  Murphy  et  al.,  2008),  and  to  the  extent  the  Rg  scatters  into  Lg, 
which  can  be  significant  (Myers,  personal  communication),  this  also  affects  the  1-2  Hz  Pg/Lg  ratios. 


Figure  6.  We  compare  1-2  Hz  and  6-8  Hz  Pg/Lg  ratios  for  NTS  explosions  averaged  from  ELK  and  KNB  to 
the  relative  CLVD  strength  determined  from  waveform  modeling  (Ford  et  al.,  2008).  Note  there 
appears  to  he  a  correlation  for  low  frequency  (excepting  METROPOLIS,  which  has  a  negative 
CLVD  relative  strength)  hut  not  for  high  frequency.  Coda-derived  spectra  show  a  peak  tied  to 
absolute  depth  (Murphy  et  al.,  2008),  which  appears  related  to  Rg  excitation  and  independent  of 
scaled  depth  of  hurial  (and  spall/CLVD).  Low-frequency  Lg  thus  appears  to  he  affected  hy  both 
absolute  depth  Rg  effects  and  relative  depth  spall/CLVD  effects. 


Amplitude  Tomography 

In  the  MDAC  corrections  discussed  so  far,  we  have  been  using  a  single  1-D  attenuation  correction  for  a  given 
station  to  a  given  region.  However,  it  is  well  understood  that  in  tectonically  complex  regions  such  as  the  Middle 
East,  there  can  be  very  significant  differences  in  the  apparent  attenuation  of  regional  phase  amplitudes  over  short 
distances.  For  example.  Figure  7  shows  a  fairly  dramatic  example  of  large  differences  in  the  relative  amplitude  of 
6-8  Hz  Pn  and  Sn  for  two  Zagros  mountain  events  of  about  the  same  size  and  distance  away.  The  top  earthquake  has 
the  large  Pn  amplitude  and  weak  Sn  amplitude  that  we  associate  with  explosions,  particularly  when  compared  with 
the  other  earthquake  (e.g..  Figure  1).  To  the  extent  these  differences  are  due  not  to  differences  in  source  depth  and 
mechanism  but  to  propagation  effects,  we  want  to  map  them  out  and  be  able  to  correct  for  them.  To  do  this,  we  are 
applying  regional  amplitude  tomography  techniques. 

Crustal  regional  phase  (Pg,  Lg)  amplitudes  corrected  for  instrument  and  an  assumed  geometrical  spreading  (A  ’)  can 
be  written  as  a  combination  of  source  (S),  path  attenuation  (R),  and  site  (P)  effects: 

-S  *R  *P 

Traditional  amplitude  tomographic  techniques  (e.g.,  Phillips  and  Stead,  2008)  can  then  solve  for  each  of  these 
terms,  with  some  tradeoffs  between  them.  However,  this  results  in  source  terms  that  can  be  significant.  Given  a  new 
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event  in  a  ealibrated  region,  we  want  to  be  able  to  apply  the  tomographieally  derived  apparent  attenuation  and  get  a 
predieted  amplitude  to  use  as  a  eorreetion  before  forming  a  diseriminant.  What  souree  term  is  assoeiated  with  the 
new  event?  One  way  to  make  the  souree  term  more  amenable  for  amplitude  eorreetion  is  to  formulate  it  in  terms  of 
the  MDAC  (Walter  and  Taylor,  2001)  modified  Brune  (1970)  souree  term: 


5,- 


FM, 

1  +  ^ 


where  co^  = 


V/3 


where  F  and  K  are  eonstants  that  depend  on  veloeity,  density,  and  average  radiation  patterns  (see  Walter  and  Taylor, 
2001,  for  details),  and  where  o  is  the  apparent  stress  in  the  region.  In  this  ease,  the  souree  term  in  the  tomography 
ean  be  reformulated  to  solve  for  the  ehange  in  moment  from  the  assumed  initial  moment,  where  the  initial  moment 
eomes  from  a  ealibrated  eoda  moment  teehnique  (e.g.,  Mayeda  et  al.,  2003)  or  other  teehnique  (e.g.,  waveform 
modeling).  In  this  ease,  for  a  new  event,  the  moment  is  first  determined  and  then  a  2D  attenuation  MDAC  eorreetion 
applied  to  get  the  predieted  amplitude  for  an  earthquake  of  that  Mw  and  loeation,  mueh  as  we  have  been  using  the 
ID  path  eorreetion  in  the  MDAC  methodology. 


Example  of  large  6-8  Hz  Pn/Sn  variation  in  the  Middle  East 


45"E  50"E  55‘E  60°E  65°E 


Figure  7.  An  example  of  large  P/S  ratio  differences  for  two  events  of  similar  size  and  distance  that  appear  to 
be  due  to  complex  path  effects.  Top  event  looks  explosion-like  relative  to  the  bottom  event.  We  are 
exploring  amplitude  tomographic  techniques  to  account  for  and  correct  2-D  path  effects  on  regional 
amplitudes  used  in  discriminants. 

Finally,  if  a  good  estimate  of  the  source  spectrum  is  available  independently,  through  a  coda  technique  or  some 
other  way,  it  is  straightforward  to  use  that  for  the  source  term  directly  and  just  solve  for  the  path  and  site  terms 
(e.g.,  Walter  et  al.,  2007b).  We  are  in  the  process  of  applying  and  testing  these  three  variations  of  amplitude 
tomography  in  the  Yellow  Sea  Korean  Peninsula  region  (see  Ford  et  al.,  these  Proceedings)  and  in  the  Middle  East. 
We  plan  to  show  some  initial  results  at  the  meeting. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Regional  discrimination  algorithms  require  calibration  at  each  seismic  station  to  be  used  for  nuclear  explosion 
monitoring.  We  apply  MDAC  procedures  to  remove  source  size  and  path  effects  from  regional  body-wave  phases. 
This  allows  comparison  of  any  new  regional  events  recorded  at  a  calibrated  station  with  all  available  reference  data 
and  models.  This  also  facilitates  the  combination  of  individual  measures  at  multiple  stations  to  form  multivariate 
discriminants  that  can  significantly  improve  performance  over  single-station  individual  measures.  We  are  working 
to  quantify  the  performance  of  P/S  ratio  discriminants  as  a  function  of  a  number  of  factors,  including  frequency, 
phase  type,  source  media  properties,  long-period  moment  tensor,  and  2-D  apparent  attenuation  patterns. 
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